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In  this  paper  a  statistical  method  for  establishing  the  minimum  energy  reservoir  needed  for  a  hybrid- 
electric  vehicle  (HEV)  is  proposed.  This  method  is  based  on  real  world  data  and  investigates  the  statistical 
properties  of  the  charge  and  discharge  events  for  a  variety  of  driving  profiles.  The  distribution  of  the  mag¬ 
nitude  of  discharge  and  charge  events  was  found  to  be  exponentially  distributed.  Using  an  exponential 
distribution  assumption,  the  probability  distribution  for  the  energy  stored  in  the  battery  was  calculated. 
This  distribution  is  a  function  of  only  two  parameters,  the  average  discharge  energy  (/xd),  and  the  ratio  of 
the  average  discharge  energy  to  the  average  charge  energy(£).  These  parameters  are  functions  of  the  drive 
profiles  of  interest  and  the  power-level  of  the  HEV  power  supply.  Based  on  this,  a  strategy  HEV  battery 
capacity  and  the  power-level  of  the  primary  power  supply  was  proposed.  This  strategy  is  of  particular 
importance  for  fuel  cell  based  HEVs  because  the  cost  of  the  fuel  cell  stack  directly  scales  with  the  power 
level. 
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1.  Introduction 

The  traditional  engineering  design  approach  is  to  design  a  sys¬ 
tem  that  meets  requirements  derived  from  a  set  of  worst  case 
scenarios.  In  design  of  hybrid  vehicles  there  is  an  additional  com¬ 
plexity  added  by  determination  of  the  degree  of  hybridization 
required  (power  ratio  between  primary  source  and  energy  reser¬ 
voir)  [1,2].  Such  designs  most  often  follow  a  “90-10  rule”— 90%  of 
the  requirements  can  be  met  by  a  design  that  requires  10%  of  the 
resources,  satisfying  the  last  10%  of  the  requirements  then  uses 
the  remaining  90%  of  the  resources.  Building  in  these  performance 
reserves,  while  necessary  in  most  cases,  adds  cost  into  the  systems. 
Therefore,  in  setting  system  requirements  it  warrants  asking  “How 
likely  is  this  limiting  requirement  occurs ?  How  much  can  we  save  if 
we  accept  a  small  probability  of  violating  this  requirement ?”  This  is 
an  approach  commonly  used  in  robust  design. 

In  this  paper  we  present  a  robust  design  methodology  that 
encompasses  statistical  aspects  of  the  system  requirements  for 
battery  sizing  in  automotive  hybrid  powertrains.  Although  not  all 
automotive  requirements  can  be  compromised,  some  performance 
attributes  do  offer  a  potential  of  large  savings  with  minor  compro¬ 
mises  in  performance.  The  method  described  here  relates  to  sizing 
of  an  energy  reservoir  which  is  depleted  through  a  random  process 
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and  cannot  be  filled  above  a  certain  level.  It  will  be  demonstrated 
through  an  example  of  sizing  hybrid  powertrains  (ICE1  or  fuel  cell 
[3-8])  for  automotive  applications. 

2.  Problem  statement,  definitions  and  assumptions 

For  this  analysis  a  simplified  hybrid  powertrain  is  assumed 
[7,9,10],  shown  in  Fig.  1,  that  consists  of  a  primary  power  source 
(internal  combustion  engine  shown  in  Fig.  1  (a)[  10]  or  fuel  cell 
shown  in  Fig.  l(b)[8]),  an  energy  storage  device  (most  likely  a 
hybrid  battery  [11])  and  a  drivetrain  (such  as  described  in  Ref. 
[12]).  The  power  needed  by  the  drivetrain  to  propel  the  vehicle  is 
either  generated  by  the  primary  source,  or  comes  from  the  energy 
storage  device,  or  some  combination  of  those  two.  Several  stud¬ 
ies  explored  this  energy  flow  as  optimization  process  with  goal 
of  improved  efficiency  and  fuel  economy  [13-15],  while  others 
focused  on  improvements  in  efficiency  and  fuel  economy  of  the 
primary  power  source  operation  [16,17].  Unlike  the  optimization 
presented  in  [  1 8  ]  no  battery  chemistry  nor  demographic  driver  data 
were  considered  here.  The  question  that  we  wish  to  answer  is :  Given 
a  distribution  of  driving  scenarios,  as  the  power  required  as  a  function 
of  time,  what  must  be  the  capacity  of  the  battery  to  meet  a  specified 
fraction  of  the  demand?  Furthermore,  we  can  ask  how  the  capac¬ 
ity  must  change  with  the  changes  in  the  power  rating  of  the  primary 
power  source? 


1  ICE:  Internal  Combustion  Engine. 


0378-7753 /$  -  see  front  matter  ©  201 1  Elsevier  B.V.  All  rights  reserved, 
doi:  1 0.1 01 6/j.jpowsour.201 1 .05.01 0 


M.  Milacic,  C.  Gearhart  /  Journal  of  Power  Sources  196  ( 201 1 )  7784-7790 


7785 


driving  wheel 


In  order  to  answer  these  questions,  several  statistical  variables 
were  introduced  (graphically  explained  in  Fig.  3): 


Tabove—  time  above  the  threshold,  P(t)  >  Ph 
Tbeiow —  time  below  the  threshold,  P(t)  <  Ph 


cabove 


^below  — 


-L 

J  lab 

L 


[P(t)  -  PH]  dt  -  energy  above PH 


P(t)dt  -  energy  below  PH 


Emi  = 


l 

J  the 


[Ph  —  P(t)]  dt  -  energy  available  to  charge 


where  P(t)  is  power  required  at  time  t,  and  Ph  is  the  maximum 
power  of  the  primary  power  source. 

Conventional  wisdom  dictates  that  the  power  rating  of  the  pri¬ 
mary  source  has  to  be  greater  than  the  maximum  power  required 
by  the  vehicle.  Questioning  this  approach  the  following  sizing  strat¬ 
egy  questions  were  posed:  How  small  can  the  primary  power  source 
be  made  without  changing  the  energy  storage  capacity  required?  If 
the  primary  source  maximum  power  is  Ph,  how  much  energy  would 
the  energy  storage  device  need  to  provide  and  for  how  long?  Can 
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Fig.  2.  Overall  drive  power  profile. 
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the  depleted  energy  from  the  energy  storage  device  (eventually)  be 
replenished? 
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3.  Analysis 

Using  data  from  the  drive  trace  shown  in  Fig.  2,  for  every  value 
of  PH ,  distributions  for  Tabove,  Pabove>  ^beiow*  and  Ffiii  as  function  of 
PH  were  generated.  These  are  shown  in  Figs.  4  and  5,  and  are  used 


Duration  of  Events 


18000 
16000 
14000 
>,  12000 
10000 
a-  8000 
ul  6000 
4000 
2000 

0 

150 


(b)  Duration  of  P<  PH 


(a)  Energy  content  while  P>  Ph 


Refill  Energy 


Energy  [kWh] 


18000 
1 6000  v 

14000 

12000 

10000. 

8000. 

6000. 

4000 

2000. 

0 

150 


100 


50 

PH  [kW] 


14 


(b)  Excess  Energy  available  while  P<  Ph 


Fig.  5.  Energy  distributions  for  the  driving  profile,  (a)  Energy  content  while  P>Ph 
and  (b)  excess  energy  available  while  P<PW. 


to  calculate  the  distribution  of  energy  associated  with  charging  and 
discharging  events.  These  distributions  are  well  approximated  by 
exponential  distributions  with  scale  parameters  /xc,  for  charging 
events,  and  /xd  for  discharging  events.  For  the  drive  traces  used  in 
this  paper,  values  of  i±c  and  /zd,  with  their  95%  confidence  intervals, 
as  function  of  Ph  are  shown  in  Fig.  6 


3.1.  Energy  transfers 

It  can  be  assumed  that  the  energy  storage  device  (battery)  is 
at  initial  state  E0  =  0.  At  any  given  time,  the  energy  storage  device 
is  either  discharging  (P>PH)  or  charging  (P<PH),  assuming  that 
threshold  level  Ph  is  the  maximum  power  that  primary  source 
can  supply.  Although  it  is  possible  to  have  times  during  which 
the  energy  does  not  change,  these  “zero  energy”  transitions  can 
be  treated  either  as  part  of  the  previous  or  part  of  the  following 
transition.  Therefore,  every  charge  event  is  followed  and  preceded 
by  a  discharge  event,  so  we  can  graphically  represent  the  energy 
history  as  shown  in  Fig.  7.  This  sequence  of  energy  level  transitions 
can  be  written  as: 


Fig.  4.  Temporal  distributions  for  the  driving  profile,  (a)  Duration  of  P>PH  and  (b)  .  .  . 

duration  of  P<PH.  Lk+1  =  (1) 
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Fig.  6.  [ic  and  iid  (with  95%  confidence  intervals)  for  various  Ph . 

where: 

f  A Ek  >  0  (charge)  k  =  2 n 
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Fig.  7.  Energy  history  with  labeled  transition  orders. 


The  odd  numbered  energy  levels  in  this  sequence  represent  the 
energy  depletion  of  the  battery  after  each  discharge  event.  The  dis¬ 
tribution  of  these  energy  levels  tell  us  how  frequently  the  battery 
is  discharged  to  a  specific  level.  To  calculate  this  distribution  (D) 
we  calculate  the  distribution  of  energy  levels  associated  with  a  sin¬ 
gle  discharge  event  (Di ).  To  this  we  add  the  distribution  of  energy 
levels  after  three  transitions  (D3)  multiplied  by  the  probability  that 
the  second  transition  did  not  replenish  the  battery  (P2)  added  by 
the  distribution  of  energy  levels  after  five  transitions  (D5)  multi¬ 
plied  by  the  probability  that  neither  the  second  (P2)  nor  the  fourth 
transitions  (P4)  replenished  the  battery,  and  so  on.  In  expanded 
form: 


D(E)  =  D,(E)  +  P2(D3(E)  +  P4(D5(E)-  •  • ))  =  j  D  2i+l  <«IM- 

1=0  \  j= 0  J 

(3) 


Now,  the  energy  level  distributions  and  energy  transition  prob¬ 
abilities  are  examined  in  more  details.  Keeping  in  mind  that  for  a 
discharge  event  AE<0,  the  energy  transition  probability  is  given 
as: 

P(AE)=  —eAE^“.  (4) 

Md 

Similarly,  the  transition  probability  for  a  charge  event  (AE>0) 
is: 

P(AE)=  —e~AE^c.  (5) 

Me 

Before  the  discharge  event  k  =  2n,  it  is  assumed  that  the  battery  is 
at  some  energy  level  E,  that  has  the  probability  distribution  D2n(E). 
After  the  discharge  event,  keeping  in  mind  that  E<  0,  the  energy 
distribution  becomes: 

o 

D2„+i(E)=  f  D2n(E)2_e(£-«)/MddE  (6) 

J  ^ 

E 


The  probability  distribution  for  E  after  the  subsequent  upwards 
transition  (step  2 n  +  2)  can  be  expressed  as: 


^2n+2(£)  =  < 


L 

/: 


D2n+^  (£)  —  e-(£-£)/^  d£  E  <  0 
Me 

D2n+-i(E)Ee-(E-E)/ncdi  £>o 
Me 


(7) 


The  corresponding  probability  that  the  subsequent  charge  will  not 
replenish  the  battery  is: 

Pm+2=  [  D2n+2(E)dE.  (8) 

J  —oo 

Normalizing  for  E  <  0  the  distribution  of  energy  levels  of  the  subse¬ 
quent  charge  is: 


D2n+2(E)=Dp+2(£). 

l2n+2 

Starting  from  (n  =  0): 


(9) 


Po=l 


and 

Di(E)=  —eE^i 
Md 
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Fig.  8.  £  =  ficIP'd  as  a  function  of  PH. 


we  can  calculate  D2  as: 

fE 

D2(E)  = 


DUE)— e(£_fi)^cd£= - ^ - ef/Pi 

Me  Me  +  Md 


and 


*-/ 


D2(£)d£  = 


Md 


Mc  +  Md 


(10) 


yielding: 

D2(£)=  —  eE/^. 

Md 

Calculating  further  D3,  D4  and  P4  we  get: 


D3(£)  = 

D4(£)  = 
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Introducing  the  parameter  £  =  Mc/Md*  shown  in  Fig.  8,  the  lowest 
order  Ps  are: 


P2  = 


P4  = 


Pe  = 


Ps  = 


Pio 

P12 


1 

f+1 

2£  +  l 

(?  +  l)2 

5£2  +4^  +  1 
(2^+  1)(4T  +  l)2 
14£3  +  14f2  +  +  1 

(f+l)2(5£2+4f+l) 

42£4  +  48£3  +  27£2  +8^  +  1 
'  (£+l)2(14£3  +  14£2  +  6f  +  l) 

132£5  +  165£4  +  110£3  +44^2  +  10£  +  1 
(f  +  1  )2(42f4  +  48f3  +  21?  +  8£  +  1 ) 


In  general,  the  probability  P2n  can  be  expressed  as  the  ratio  of 
two  polynomials  of  f : 


P2n 


•^2  n(£) 
?2n(?)’ 


(11) 


For  a  specific  hybrid  vehicle  configuration,  the  parameter  £  is  a 
function  of  P#,  as  shown  in  Fig.  8  for  our  example  vehicle.  This 
parameter  provides  insight  into  the  energy  balance  between  dis¬ 
charge  and  charge  events;  higher  values  of  f  means  that  the  energy 
used  during  discharge  is  more  likely  to  be  replenished  in  subse¬ 
quent  charge  event.  In  contrast,  smaller  values  of  £  (such  as  f  <  1) 
means  that  energy  discharged  is  less  likely  be  replenished.  If  £  is 
too  small,  the  charging  events  will  not  be  able  to  keep  up  with  the 
discharging  events  and  the  battery  will  eventually  become  totally 
discharged.  This  will  occur  if  the  primary  source  power  (low  Ph) 
rating  is  too  low  for  the  drive  traces  of  interest. 

If  parameter  a2n  is  used  to  represent  the  coefficient  of  £n_1  in 
the  polynomial  J\f2n ,  then  the  inductive  formula  for  P2n  that  is  true 
for  all  values  of  n  >  1  can  be  written  as: 


P2n+2  =  1  - 


C*2n£n+1 

(?  +  l  )++,(  + 


(12) 
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The  parameter  a2n  is  discussed  further  in  Section  4. 
The  first  several  orders  of  D2n+1  are: 

eE/ Md 
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The  distributions  D2„+i  can  be  calculated  (n  >  1 )  using: 


f^2n+l(£)  — 


(-i)nE0^d  r  k 


(n  +  k-  1)! 
k!(n  —  /<  —  !)! 


(E(f + i))n-fcV5ffe 


(13) 


The  first  few  P  parameters  as  a  function  of  Ph  calculated  based 
on  drive  profile  for  our  example  vehicle  are  shown  in  Fig.  9. 

Going  back  to  Eqs.  (3 )  and  (12),  and  realizing  that  after  canceling 
terms  from  the  multiplication  of  the  P  terms 


nw 

j= i 


•A/2n(£)  1 

(f+l^-^ff+l)"' 


(14) 


where N2n{$)  is  defined  in  Eq.  ( 1 1 ).  We  note  that  D(£ )  can  be  written 
as: 


D(E)  =  D1(£)  +  ^ 

n= 1 


(^D2n+1(E) 


J^2niX)  \ 
(4:  + 1)2"-1  /  ■ 
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P2n  =  f(PH> 


0.68269  after  68.269%  of  discharge  events.  We  can  see  that  in  this 
case,  when  the  primary  power  supply  is  larger  than  75  kW  the  con¬ 
tours  are  nearly  vertical.  This  means  that  the  size  of  the  primary 
power  supply  has  almost  no  impact  on  the  distribution  of  energy 
in  the  battery. 

It  can  also  be  seen  that  for  PH^  lOkW  (f  =  1  point)  the  contour 
lines  become  horizontal  meaning  that  the  energy  capacity  require¬ 
ment  increases  (almost)  indefinitely.  As  mentioned  before,  f  =  1  is 
limiting  case  below  which  the  primary  power  source  is  not  capable 
of  replenishing  the  energy  used. 

4.  Parameter  a 

Introduced  in  Eq.  (12)  as  the  coefficient  associated  with  n-  1 
power  of  £  in  the  polynomial  A/E.  the  parameter  ot2 n  deserves  a 
closer  look,  especially  since  it  is  not  dependent  on  £.  From  Eq.  (12) 
for  n  >  1 : 

T>2n+2(t)  =  (t+P2^2  n(4T) 

and 


Fig.  9.  P2n  as  a  function  of  PH. 


Using  Eq.  (13)  and  rearranging  the  terms  it  becomes: 


D(E)  = 


eE/iid 

Md 


n-1 

k= 0 


(n  +  fc-p!  (-?)* 
k\n\(n-k-  1)!  (£  +  i)n+k 


After  normalization,  the  normalized  distribution  D(E)  can  be 
written  as: 


D(E)  = 


eE/tid 


oo 


n-1 
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<-»"E  [ 
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M2n(Q 
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The  cumulative  distributions  of  D(E)  as  function  of  Ph  for  our 
demonstration  vehicle  are  shown  in  Fig.  10.  Each  contour  on  this 
graph  represents  a  specific  probability  that  the  battery  will  reach  a 
certain  energy  level  after  a  discharge  event.  For  example  the  energy 
level  in  the  battery  will  be  below  or  to  the  right  of  the  line  marked 


Probablility  Contours 


■A/2n+2(0  =  ^2n+2(0  “  <*2 n£n+1  • 

Therefore,  a  recursive  calculation  for  A/E(f).  with  initial  values  for 
n  =  0,  A/E?)  =  1,  and  a0  =  1,  can  be  expressed  as: 

AWf)  =  (?  +  1  )2A /2„(?)  -  a2nr+1  •  (16) 

Now,  denoting  the  coefficients  of  the  polynomial  A/E(?)  as  ani-  with 
i  =  0. ..n-1, 

n-1 

•^E(?)  =  y^an,n-k-l?k- 

k= 0 

Noting  that  a2n  =  an,o ,  the  dependency  of  the  coefficients  of  the 
A/*2n+2 ( 0  can  be  written  as: 

t*n+i,j  —  tzn+i  j_i  +  2anj  +  onj+ 1 

with  the  exception  of  (after  term  o(2n?n+1  is  canceled): 
fln+ 1,0  =  C^2n+2  =  2an?0  +  an,l  =  2  a2n  +  On,l-  (l7) 

Defining  the  tridiagonal  n  x  n  matrix  0  (0  =  tridiag(l,  2,  l))as: 
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the  coefficients  of  polynomial  A/E(0  can  be  calculated  as: 


<3n,0 

"1" 

Gn,l 

0 

-<3n,n-l  - 

_0_ 

Calculating  powers  of  0  using  methods  described  in  Refs.  [19,20], 
is  straight  forward.  Of  particular  interest  is  a2n  =  an  0  that  can  be 
expressed  as: 


4(2n-l)  4n+1(2n  -  1)! ! 

“2n  _  2n  +  2  “2n_2  _  2(2n  +  2)! ! 


(19) 


Some  coefficients  that  are  easily  calculated  are  an<n_ i  =  1  and 
Q-n,n- 2  =  2  n. 


Fig.  10.  D(E)  for  various  Ph. 
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5.  Special  case  if  =  1 


This  case  is  singled  out  because  it  is  the  limiting  case  between 
converging  and  diverging  energy  distributions.  For  this  value  of  £ 
the  distributions  for  charge  and  discharge  are  the  same.  In  this  case 
the  energy  state  may  eventually  return  to  its  initial  state,  but  the 
number  of  high  order  excursions  becomes  significant,  and  mathe¬ 
matically  speaking  convergence  of  D(E)  cannot  be  guaranteed.  For 
that  reason,  it  deserves  closer  look. 

Using  f  =  1,  we  can  expand  the  transition  probabilities  P  as: 


P;-2 


126 


_462 

10  140  12  504 


or  after  cancelation  (in  each  case  by  the  previous  a2n  parameter,  as 
described  in  Eq.  (12))  the  P’s  become: 


P2  =  2 
Pw  =  W 


which  can  be  written  as: 

P2n  =  1  -  i  (20) 


An  inductive  proof  is  straight  forward,  using  Eqs.  (12)  and  (19). 
Using  Eq.  (12)  we  can  write: 


n  — 


4A/i„(l) 
(2n  +  2) 


(21) 


and  from  Eqs.  (19)  and  (21)  we  get: 


A/2n+2(l)  = 


4n 


(2n  +  l)M 
(2n  +  2)M 


(22) 


6.  Summary  and  conclusions 


native)  and  battery  capacity  based  on  an  acceptable  reliability  limit. 
For  example,  the  thick  line  in  Fig.  10  represents  a  reliability  level  of 
only  3  instances  out  of  one  million  discharge  events  would  exceed 
the  battery  capacity.  In  an  example  of  such  a  scenario,  instead  of 
providing  120kW  for  1  s  the  drivetrain  would  receive  80  kW  for 
1.5  s,  which  can  be  hardly  noticeable  by  the  driver. 

Finally,  it  should  be  pointed  out  again  that  this  proposed 
sizing  method  for  hybrid  powertrain  is  purely  mathemati¬ 
cal  and  does  not  take  into  consideration  neither  physical 
limitations  (chemistry,  capacitance,  losses,  degradation  over 
time,  etc.),  nor  vehicle  requirements  (acceleration,  gradability, 
etc.). 
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Based  on  Fig.  10,  it  is  possible  to  approximate  the  minimum 
engine  power  (conventional  internal  combustion  or  any  other  alter- 


